Double Exchange Ferromagnetism in the Peierls Insulator State 



S. Nishimoto 1 and Y. Ohta 2 
1 Institut fiir Theoretische Festkorperphysik, IFW Dresden, 01171 Dresden, Germany and 
2 Department of Physics, Chiba University, Chiba 263-8522, Japan 
(Dated: July 9, 2012) 

We study the effects of opening of the band gap on the double exchange ferromagnetism. Applying 
the density-matrix renormalization group method and an analytical expansion from the dimer limit 
to the one-dimensional double exchange model, we demonstrate for a relevant region of the exchange 
coupling that, in the weak dimerization regime, the Peierls gap opens in the fully spin-polarized con- 
duction band without affecting its ferromagnetism, whereas in the strong dimerization regime, the 
ferromagnetism is destroyed and the Mott gap opens instead, leading the system to the antiferro- 
magnetic quasi-long-range order. An insulator version of the double exchange ferromagnetism is 
thus established. 
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Magnetism and electronic transport properties of ma- 
terials are closely related to each other, e.g., insulat- 
ing transition-metal oxides are typically antifcrromag- 
netic and ferromagnetism usually goes hand in hand with 
metallicity [l[ , and elucidation of the mechanisms of this 
relationship is one of the major issues in the field of 
physics of strong electron correlations. A well-known ex- 
ample is the double exchange ferromagnetism that occurs 
in mixed systems of localized spins and itinerant elec- 
trons interacting via the Hund's rule coupling, where the 
coherent motion of the itinerant electrons aligns the lo- 
calized spins ferromagneticalby to gain in kinetic energy 
of the itinerant electrons [2|-|4[ . 

A question then arises as to what happens in this fer- 
romagnetism if the coherent motion of electrons ceases 
due, e.g., to the opening of the band gap. This issue 
comes into a real question in ferromagnetic hollandite 
K 2 Cr 8 0i6 @, where the double exchange mechanism in- 
duces the three-dimensional (3D) full spin polarization 
in the system Q , and then the metal-insulator transition 
follows in its fully spin-polarized quasi-one-dimensional 
(ID) conduction band by the Peierls mechanism 0,11], 
without affecting its 3D ferromagnetism. Thus, the un- 
common ferromagnetic insulating (FI) state is realized in 
this material 0. 

A naive answer to the above frequently asked question 
may then be that the ferromagnetism can survive if the 
band gap is small enough in comparison with the width of 
the conduction band and therefore the motion of conduc- 
tion electrons, though not coherent, is not significantly 
suppressed. However, to the best of our knowledge, no 
quantitative theoretical studies have been made on this 
issue whether this is actually the case. 

In this paper, we address this issue from the theoret- 
ical standpoint. We apply the numerical density-matrix 
renormalization group (DMRG) technique @ and ana- 
lytical expansion from the dimer limit to the ID dou- 
ble exchange model that simulates the quasi-lD chain of 
K2CrgOi6 and study the effects of opening of the band 



gap on the double exchange ferromagnetism. We calcu- 
late the total-spin quantum number and charge gap of 
the system and extract the ground-state phase diagram 
of the model. 

We will thereby demonstrate for a relevant region of 
the exchange coupling that, in the weak dimerization 
regime, the Peierls gap opens in the fully spin-polarized 
conduction band without affecting its ferromagnetism, 
whereas in the strong dimerization regime, the ferromag- 
netism is destroyed and the Mott gap opens instead, due 
to the effective "on-dimer" Coulomb interaction, which 
leads the system to the antiferromagnetic quasi-long- 
range order. The metallicity itself is not therefore a 
necessary condition for the realization of the double ex- 
change ferromagnetism and thus an "insulator version" of 
the double exchange ferromagnetism is established. This 
is a route to realization of insulating fcrromagnets, dif- 
ferent from those discussed in doped LaMnOj where the 
orbital ordering plays an essential role 13- I1J ■ 

Our model (see Fig. [IJ contains the terms of Peierls 
dimerization and Hund's rule coupling and is defined by 
the Hamiltonian 
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where c\ a is the creation operator of an electron with 
spin a (=t,4) at site i, Si is the spin operator of a 
conduction electron at site i, t is the vector of Pauli 
matrices, and Si is the quantum spin operator (of spin 
1/2) of a localized electron at site i. The hopping pa- 
rameter between the nearest-neighbor sites is defined as 
ti,i+i = [l — ( — l) l 5/2]t with the dimerization strength 
5 of < 8 < 2; we in particular define t± = (1 ± S/2)t. 
Jh is the strength of the Hund's rule coupling. L is the 
number of sites in the system, where the site contains a 
conduction orbital and a localized spin. We confine our- 
selves to the case at quarter filling of conduction electrons 
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FIG. 1: (Color online) Schematic representations of the dou- 
ble exchange model in (a) the weak dimerization regime and 
(b) strong dimerization regime. In (c), we illustrate the non- 
interacting band structure e(k) of our model with the lattice 
dimerization S in the unfolded Brillouin zone. 



n = N/L = 1/2, where N is the number of conduction 
electrons in the system. We introduce 6 as a controllable 
input parameter rather than the order parameter of the 
spontaneous lattice dimerization since our purpose is to 
study the effects of 8 on the electronic states of the model. 

We use the DMRG method to investigate the ground- 
state properties of the system Eq. (fl]) , where the open- 
end boundary conditions are applied. We study the 
model with several lengths of L = 8 to 20 with keeping 
to = 1200 to 3200 density-matrix eigenstates in the renor- 
malization procedure; in this way, the largest truncation 
error, or the discarded weight, is of the order of 10 -11 . 
Note that we have to keep relatively-numerous density- 
matrix eigenstates to extract the true ground state from 
a number of nearly degenerate magnetic states, in par- 
ticular, in the vicinity of the phase boundaries. The ex- 
trapolation to the thermodynamic limit L — > oo is made 
in the results presented unless otherwise indicated. 

First, let us present the total-spin quantum number 
S of the system, which is calculated directly from the 
ground-state expectation value of S 2 defined as (S 2 ) = 
S{S+1) = Zij ({svS j }+2{s r S j } + {SrS j }). The results 
are given in Fig. [2] as a ground-state phase diagram. We 
find that there are three phases, S = 0, < S < S max , 
and S = S max , depending on the values of Jh and S, 
where S max denotes the full spin polarization. As ex- 
pected, and in agreement with the previous calculations 
at S = HI- HI], the S = S max phase appears when 
Jh is large and S = phase appears when Jh is small. 



FIG. 2: (Color online) Calculated ground-state phase dia- 
gram of the ID double exchange model at quarter filling with 
the lattice dimerization S. We find the FI (ferromagnetic in- 
sulating), PI (paramagnetic insulating), FM (ferromagnetic 
metallic), and PM (paramagnetic metallic) phases. The FI- 
PI phase boundary at S — > 2 is determined by the analytical 
expansion from the dimer limit and is shown by the open di- 
amond and dashed line. The FM and PM phases at S — or 
Jh = are also indicated by the thick lines. 



polarization also appears as in Ref. [15| at S = (see also 
Fig. EJa) below). These phases are retained even when 
the dimerization 5 is introduced. We note that the two 
critical values of Jh that separate between the S = S roax 
and S = phases increase with increasing S and that the 
region with the intermediate spin polarization becomes 
narrower and vanishes at 6 — !> 2 (see below). 

Next, let us calculate the charge gap A defined as A = 
lim^oo A(L) with A(L) = \[E^ +2 {L) + E^ 2 (L) - 
2Eq(L)] , where Eq(L) is the ground-state energy of the 
system of size L with N electrons. The gap A is defined 
with the prefactor 1/2, so that the single-particle band 
gap is identical to A in the present case where pairing 
interactions are absent. The results are shown in Figs. [3] 
(b)-(d). We find that the charge gap opens in the entire 
parameter space except at the lines 5 = and Jh = 0. 
The model at Jh = is trivial, where the conduction 
electrons, decoupled completely from the localized spins, 
behave just as the noninteracting electrons, resulting in 
the PM phase. 

The model Eq. (Q]) at 6 = on the other hand is highly 
nontrivial and much work has been done in recent years 
[UEBl : basically, there appears the ferromagnetic metal- 
lic (FM) phase for the large Jh region, which changes into 
the paramagnetic metallic (PM) phase when Jh becomes 
small. In addition, it has been claimed jljl 113 that 
And in between, the phase with the intermediate spin there appears the region of phase separation in particu- 
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lar near half filling and the "spiral" phase with a long- 
wavelength antifcrromagnctic correlations at the FM-PM 
phase boundary. Although we have not detected any in- 
dications of the phase separation at least at quarter fill- 
ing in our accurate DMRG calculations with very large 
m values, our results obtained are consistent with the 
results of the previous work EMU: we find that cither 
the FM phase (when Jh is large) or PM phase (when Jh 
is small) is realized, and in between there is the partially 
spin-polarized metallic phase (see Figs. [2]and[3][a)) that 
may correspond to the spiral phase predicted in Ref. [TBI . 




FIG. 3: (Color online) (a) Calculated normalized total-spin 
quantum number S/Smax as a function of Ju/t- Also shown 
are the calculated results for the charge gap A: (b) <5 depen- 
dence in the FI phase, (c) 8 dependence in the PI phase (where 
the dotted line indicates the result of the strong-dimerization 
expansion), and (d) Jn/t dependence in the PI phase. Inset 
of (b) and (c) shows examples of the finite-size scaling of the 
gap. In (e) and (f), the calculated spin structure factors S(q) 
in the PI phase without and with the lattice dimerization, 
respectively, are shown. 



Let us then introduce the lattice dimerization S > 0. 
The results are the following: In the S = S max region, 
we find that the charge gap of A/t = 26 opens as shown 



in Fig. Gl^b). This can simply be understood because in 
this region we have the noninteracting band of spinlcss 
fermions at half filling, and therefore the lattice dimer- 
ization 5 opens the band gap of the size A/t = 26. The 
FI phase is thus realized. Since the FM phase due to 
the double exchange mechanism is continuous to this FI 
phase, we may naturally refer to it as the insulator ver- 
sion of the double exchange ferromagnctism. 

In the S = region, we find that the gap actually 
opens as A cx 6 in the weak dimerization limit as shown 
in Fig. EJc). The size of the gap increases as the value 
of Jh increases as shown in Fig. [3Jd) . This phase with 
A > may then be denoted as the paramagnetic insu- 
lating (PI) phase. The Fourier transform of the spin- 
spin correlation function for the conduction electrons 
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Si ■ Sj) (as well as that for the 
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localized electrons, see Fig. 2 of Refill 5l) may characterize 
this phase. The calculated results are shown in Figs. [3] 
(e) and (f), where we find that the antifcrromagnctic spin 
correlation of the wavevector of q = it/2 is enhanced and 
that the lattice dimerization inducing the localization of 
conduction electrons further enhances this correlation. 
This results may therefore be interpreted as an enhance- 
ment by the lattice dimerization of the "island" state 
predicted in Ref. [H, where the high-spin (S = 3/2) clus- 
ters, formed by a conduction electron coupled ferromag- 
netically with the two neighboring localized spins, are 
arranged antiferromagnetically with the quasi-long-range 
order [Til, E3- I n higher spatial dimensions, this phase 
may well fall into the true long-range antifcrromagnctic 
order, resulting in the antifcrromagnctic insulating (AFI) 
phase. This situation resembles that of the "dimer-Mott" 
phase [i(J [HJ in the dimerized Hubbard model at quarter 
filling although in the latter the spins are of S = 1/2. 

In the < S < S max region, the charge gap also opens 
at 5 > 0, where its size increases rapidly with increasing 
Jn/i. We thus have the FI phase here as well, which may 
be the insulating spir al p hase continuous to the metallic 



one predicted in Ref. |15 . 

Now, let us discuss the strong dimerization limit, 
where we start with the highly correlated clusters C; 
(I = 1, • • ■ ,L/2) coupled weakly to each other through 
the hopping parameter t_ (see Fig. [ljb)). Each of the 
clusters consists of the two conduction orbitals and two 
localized spins. In the ground state, the cluster contains 
one conduction electron (and two localized spins) and the 
internal three spins are fully polarized. The lowest en- 
ergy of the single cluster is e(l) = — Jh/4 — t + , where 
the conduction electron is in the bonding state of the 
two conduction orbitals. The eigenstates of the cluster 
with the and 2 conduction electrons are also derived 
exactly and the lowest energies are found to be e(0) = 

and e(2) = — \J + 16i+/2, respectively. In the strong- 
dimerization limit of the PI phase, we can therefore map 
our system Eq. ([T]) onto an effective single-band Hub- 
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bard model defined in terms of the bonding orbital on 
each dimcr. The Hamiltonian may be written as 

L/2 L/2 

H cff = t eff ( b t* b i+h° + H - c + UcS 4t n U> ( 2 ) 

i=l,<r i=l 

with the creation operator of an electron on the bond- 



ing orbital b\ = (c|j 



„t 



a )/y/2 and 



6t 6,; 



We obtain the effective hopping integral i c rr = t-/2 
and effective "on-dimer"" Coulomb interaction U c s = 
e(2) + e(0) — 2e(l). In this mapping, the localized spins 
contribute only to f/ e ff and their degrees of freedom are 
not explicitly involved in the operator b\ . An analytical 
expression for the charge gap may thus be derived as 

A = U e «-2t- 



= - y/jg + 4i 2 (2 + (5) 2 /2 + 2t8 + J H /2, (3) 

which is valid up to the first order of 2—8. The result thus 
obtained is shown in Fig. (3{c) as a dotted line, where we 
find that the agreement with the DMRG result is very 
good. In the FI phase, on the other hand, the charge 
gap is always A/t = 25 in its entire region, independent 
of Jh- The two gaps are therefore discontinuous at the 
phase boundary. 

Next, we derive the effective exchange interaction J e g 
between conduction electrons on the neighboring clus- 
ters. Wc first take a direct product of the isolated clus- 
tcrs (Ilj=i Ci) t° be the unperturbed ground state of the 
system. Then, taking into account all the processes up 
to the second order of the hopping i_, we obtain the 
expression for J e g as 
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where we should note that, if only the lowest intermedi- 
ate state is taken into account in the second-order pro- 
cess, we obtain J cff = 4i 2 _/(7 of f, which is always positive 
(or antifcrromagnctic) . Wc thus find that, depending on 
the value of J e rr (either positive or negative) , the ground 
state realized is either the PI phase or the FI phase. The 
critical value of Jh at 8 — > 2 is found to be 29.004i, which 
determines the FI-PI phase boundary at 8 — > 2. The re- 
gion of the intermediate spin state does not appear here. 
The result for the phase boundary obtained from Eq. ^ 
is shown as a dashed line in Fig. [21 where we find again 
that the agreement with our DMRG result is very good, 
reinforcing the validity of our phase diagram of Fig. [2j 

In summary, we studied the effects of opening of the 
band gap on the double exchange ferromagnetism. Wc 
applied the DMRG technique and analytical expansion 
from the dimcr limit to the ID double exchange model 



at quarter filling with lattice dimcrization and obtained 
the ground-state phase diagram. We found three phases: 
the FI phase with the Peierls gap and full spin polariza- 
tion, the PI phase with the Mott gap and dominant an- 
tiferromagnetic spin correlations, and the FI phase with 
the partial spin polarization. The results for Jn/t > 15 
demonstrated that, in the weak dimcrization regime, the 
Peierls gap opens in the fully spin-polarized conduction 
band without affecting its ferromagnetism. Therefore, 
the metallicity itself is not a necessary condition for 
the realization of the double exchange ferromagnetism. 
The concept of the insulator version of the double ex- 
change ferromagnetism was thus established. In the 
strong dimcrization regime, on the other hand, the ferro- 
magnetism is destroyed at Jn/t < 29 and the Mott gap 
due to the effective on-dimer Coulomb interaction opens 
there with the antiferromagnetic quasi-long-range order 
in the system. 

The uncommon FI state realized in E^CrgO^ 0] then 
means that this material is in the weak dimcrization 
regime with the Peierls gap and full spin polarization. 
A recent experiment suggests [13] that, by applying high 
pressures of > 2 GPa, the FM phase is suppressed very 
rapidly, while the metal-insulator transition remains al- 
most unchanged, leading to the transition from the PM 
phase to the PI or AFI phase by lowering temperature. 
It may then be quite interesting to point out that, if the 
applied pressure decreases the value of Jn/t, this might 
correspond to the phase change in the quasi-lD chains 
from the FI to PI phase as in our phase diagram given 
in Fig. [2j where the intermediate spin state, or the spiral 
state of Ref. Il5l . may also be predicted to appear under 
high pressure. We hope that our work presented here will 
stimulate further searches for new phenomena and ma- 
terials with intriguing magnetic and transport properties 
derived from the interplay between the double exchange 
and Pcierls/Mott mechanisms. 
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